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ABSTRACT 

We identify a new, nearby (0.5 < d < 10 kpc) stream in data from the RAdial Velocity Experiment 
(RAVE). As the majority of stars in the stream lie in the constellation of Aquarius we name it the 
Aquarius Stream. We identify 15 members of the stream lying between 30° < I < 75° and —70° < 
6 < —50°, with heliocentric line-of-sight velocities V|os ^ — 200kms~^. The members are outliers 
in the radial velocity distribution, and the overdensity is statistically significant when compared to 
mock samples created with both the Besangon Galaxy model and newly- developed code Galaxia. The 
metallicity distribution function and isochrone fit in the log g - Tos plane suggest the stream consists of 
a 10 Gyr old population with [M/H] ^ —1.0. We explore relations to other streams and substructures, 
finding the stream cannot be identified with known structures: it is a new, nearby substructure in 
the Galaxy's halo. Using a simple dynamical model of a dissolving satellite galaxy we account for 
the localization of the stream. We find that the stream is dynamically young and therefore likely 
the debris of a recently disrupted dwarf galaxy or globular cluster. The Aquarius stream is thus a 
specimen of ongoing hierarchical Galaxy formation, rare for being right in the solar suburb. 
Subject headings: Galaxy: halo - Galaxy: kinematics and dynamics - Galaxy: solar neighbourhood 



1. INTRODUCTION 

Under the current paradigm of galaxy formation galax- 
ies build via a hierarchical process and our Galaxy is 
deemed no exception. Relics of formation are observed 
as spatial and kinematic substructures in the Galaxy's 
stellar halo. Recent observations such as those from the 
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Sloan Digital Sky Survey (SDSS) have brought a large in- 
crease in the detections of substructures within the outer 
reaches of the halo (out to d < 80 kpc) . These streams 
have usually bee n detected as spatial overdcnsitie s from 
photometry e.g., Yanny et al.| ( 2000[); [Mai cws ki et al. 
( 2003| ; IBeiokurov et al.| ( |:j006p ; INewbcr g et~a:L| ^' pUg^ 
Ivlany of these structures have been identified as be- 
longing to the debris of the Sagittarius dwarf spheroidal 
galaxy (Sgr dSph), which traces the polar orbit of this 
galaxy as it merges with the Milky Way. Fur thermore, 
after subtracting such prominent substructures [B"ell et al.| 
(2008) observed a dominant fraction of the halo to de- 
viate from a smooth distribution, consistent with being 
primarily accretion debris. 

Closer to the Sun the spatial coherence of streams 
and substructures is not so easily discernible and most 
streams of stars are visib le only as velocity stru ctures, 
such as the Helmi et al. ( )1999[ ) stream. Indeed, Helmi 
(20091 has shown that only at distances greater than 
~ 10 kpc do we expect that the structures associated 
with tidal debris to be observable as spatial overdcnsi- 
ties. Therefore, if we wish to identify and study struc- 
tures within the inner reaches of the halo - where they are 
most accessible for high resolution follow-up observations 
- we must search utilizing kinematic data. 

Kinematic surveys of the solar neighbourhood are 
therefore ideal to detect substructures in the nearby re- 
gions of the Galaxy's halo. RAVE (RAdial Velocity Ex- 
periment) is an ambitious program to conduct a 17,000 
square degree survey measuring line-of-sight velocities, 
stellar parameters, m etallicities and abunda nce ratios of 
up to 1 million stars ( Steinmetz et al.||2006 ). RAVE uti- 
lizes the wide field (30 deg^) multi-object spectrograph 
6dF instrument on the 1.2-m UK Schmidt Telescope of 
the Anglo- Austrahan Observatory (AAO). RAVE's in- 
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put catalogue for the most parl[^ has only a magnitude 
selection criterion of 9 < / < 13, thus creating a sample 
with no kinematic biases. The observations are in the Ca- 
triplet spectral region at 840 nm to 875 nm with an effec- 
tive resolution oi R ~ 7500. Starting in April 2003, at the 
end of 2009 RAVE had collected more than 400,000 spec- 
tra, rave's radial velocities are accurate to 1.3 km s^^ 
when compared to external measurements, while t he re- 
peat obse rvations exhibit an accuracy of 2kms^^ (jZwit-J 



ter et al. 2008). These highly accurate radial velocities 
make rAvE ideal to search for kinematic substructures 
in an extended region around the sun. Indeed, with 
RAVE we now move away from s t udyin g the solar neigh- 
bourhood (e.g. Nordstrom et al. (2004): d < 0.2 kpc) to 
examining the solar suburb [d < 4 kpc) . 

Using rave's highly accurate radial velocities, we 
have discovered a stream that lies mostly within the con- 
stellation of Aquarius at a distance of 0.5 < d < 10 kpc, 
in the direction {I, b) ^ (55°, —60°) and at Vios = 
— 200kms~^. The velocity places the stream as part of 
the Galaxy's halo. As it lies in the direction of the con- 
stellation of Aquarius we have named it the Aquarius 
stream. The detection of this stream is described in Sec- 
tion [2] In Section [3] we compare the RAVE data to mock 
data from the Besangon Galaxy mode and the newly- 
developed galaxy modelling code Galaxia, which offers 
a number of significant advantages. Using these models 
we determine the significance of the detection and con- 
strain its localization. In Section|4]we use RAVE's stellar 
parameters combined with 2MASS [JHK] photometry 
to infer basic properties of the stream population and 
derive distance estimates. We also use Reduced Proper 
Motions to obtain another estimate of the distances. The 
stream appears to be highly localized on the sky which 
is interesting considering the apparent proximity of the 
stream. In Section [5] we explore possible connections of 
the Aquarius stream to other known spatial and kine- 
matic streams, finding that it is not linked to any pre- 
viously reported structure. In Section [6] we investigate 
possible connections to other (marginal) over-densities in 
the RAVE dataset, and conclude that the stream is un- 
likely to be associated with any of them. A simple model 
of the recent disruption of a satellite in the Galaxy's po- 
tential is able to account for the observed localization. 
The Aquarius stream thus is a new and nearby enigma 
in the Milky Way's halo. 

2. DETECTION IN RAVE 

2.1. The sample 

RAVE measures the velocities of stars that are selected 
purely on the basis of their photometry, so it is free of 
kinematic biases. Over most of the sky the probability of 
a star's selection depends entirely on its apparent mag- 
nitude; only in directions towards the Galactic Centre is 
selection based on colo ur as w ell as magnitud e (DRl: 
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Steinmetz et al. ( |2006| ), DR2: _ 

I'urthermore, rAVE's radial velocities are accurate to 
< 2kms~^ so fine substructures are best detected using 
radial velocities alone: combining them with proper mo- 

Red giants in the direction of rotation were also targeted 
between 225° < I < 315°, 5° < |b| < 25° with J K > 0.5. This 
region is not discussed in this paper however. 



tions and distances mean a significant loss of accuracy. 
The Aquarius stream was discovered in RAVE data as 
a structure seen in heliocentric radial velocity vs Galac- 
tic latitude/longitude space. When the stream was first 
noted, it was found to be most clearly defined by faint 
stars with low gravities, which suggests that the structure 
is at some distance from the Sun. Removing foreground 
giants enhances its visibility. 

We use the internal release of RAVE from January 2010 
that contains 332,747 RVs of 252,790 individual stars. 
We use only those observations for which the signal-to- 
noise ratio SNR > 13 and the Tonry and Davis cross- 
correlation coefficient i? > 5 to remove potentially er- 
roneous observations. Note that, since not all observa- 
tions have the more accurate signal-to-noise estimation, 
5*2^^, we use the SNR value which can underestimate 
the signal-to- noise (see DR2). For multiple observations 
of single stars the T^os were averaged, as were the stellar 
parameters for those observations that yielded an esti- 
mate of these parameters. 

The Aquarius stream was found in the Galactic lati- 
tude slice —70° < 6 < —50°. As described above, it is 
also more marked for fainter stars. We therefore intro- 
duce an upper brightness limit to enhance the visibility 
of the stream. As noted in the first and second data re- 
lease papers, a subset of the RAVE input catalog have I 
magnitudes from the SuperCOSMOS Sky Survey (Ham- 
jblyi'20011) , which show an offset to DENIS / magnitudes. 
Not all RAVE stars have DENIS / magnitudes either. 
We therefore turned to 2MASS bands for our magnitude 
limit, even though this tends to bias against cool stars 
in our sample, and we potentially miss some candidates. 
We found that a Hmit of J > 10.3 produced the best 
differentiation of the stream from the background popu- 
lation, removing the brighter, nearby giants. 

2.2. Detected overdensity 

Figure [1^ shows the structure seen in heliocentric ra- 
dial velocity, Vios, against Galactic longitude, I, for the 
stars with the selection criteria —70° < b < —50°, 
J > 10.3. A clear structure begins at Vios ^ — 150kms~^ 
at I = 30° and extends down to T^os = — 200kms~^ at 
I — 75°. This overdensity is particularly clear in Figure 
[T]d, where we plot the histogram for V\os in the region 
-70° < 6 < -50°, J > 10.3, 30° < / < 75°. The stream 
can be seen as an excess of stars at negative velocities 
that is distinct from the general population. 

We establish limits of -250 < Vios < -150kms^\ 
30° < ? < 75°, J > 10.3 to choose 15 candidates of 
the Aquarius stream, which are outlined by the red box 
in Figure [l] and listed in Table [l] Many stream candi- 
dates lack stellar parameter estimates, since they were 
observed early on by RAVE (DRl does not include such 
estimates; see the data release papers for details). The 
average SNR is 20 for the stream candidates and 1 star 
(C2234420-082649) has a repeat observation, which is 
listed to show the consistency of the Vios results. As 
a double-check, the template fits to each of the spectra 
were eye-balled as were the zero-point fits (using sky ra- 
dial velocities) for the fields the stars were observed in. 
No abnormalities were detected. 

The RAVE internal release includes PPMX proper mo- 
tions (Roeser et al. 2008, ). However, for our stream 
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Figure 1. (a) Vio^ as a function of galactic latitude for RAVE data with —70 < b < —50, J > 10.3. The Aquarius Stream is identified as 
an overdensity of stars with —250 < Mos < — 150kms~^, 30° < I < 75°, as delimited by the red box. (b) The histogram of Vios with the 
additional constraint 30° < I < 75° clearly shows the stream as an anomalous feature in the wings of the velocity distribution. The grey 
shading displays the itlcr limits. 
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Table 1 













The Aquarius Stream candidates selected from the RAVE data and their parameters. The proper motions are from PPMXL. 



candidates we u se in the following analysis PPMXL 
proper motions (Roeser et al. 2010)), where the aver- 
age proper motion error tor the stream stars is reduced 
from = 6.8masyr~^ in PPMX to = 4.3masyr~^ in 
PPMXL. These proper motions are also listed in Table 

m 

The average heliocentric radial velocity of the stream 
is V\os — — 199± 27knis^^ and its Galactocentric radial 
velocity, i.e. the line-of-sigh t velocity in the Ga lactic rest 
frame (see equation 10-8 of Binney & Tremaine (1998 )), 
is V„ 



gal 



-93 ± 25kms . When compared to Vio 



-120 ± 100kms"\ Fgai = ± 100 kms"^ for the halo 

and Vios = -30±45kms~\ Fgai = 90±45kms"^ for the 
thick disk at (1, b)=(55°, —60°), this velocity indicates 
that the group to be a halo feature. However, it still has 
quite a large velocity even for the halo. 

3. MODEL COMPARISONS 

3.1. Besangon and Galaxia models 

To establish the statistical significance of the Aquarius 
overdensity we compare the RAVE sample to moc k sam 
pies created using the Besangon Galaxy model (Robin 
et al. 20031 a nd the ne wly developed galaxy modeling 
code Galaxia ( Sharma et al. 2010i ). Galaxia is based on 



the Besangon Galaxy model, but with several improve- 
ments. The first is a continuous distribution created 
across the sky instead of discrete sample points. Sec- 
ond is the ability to create samples over a n angular area 
of arb itrary size. Third, it utilizes Padova ('Girardi et al. 



2002 1 isochrones which offer support for multiple photo 
metric bands. Fourth, Galaxia offers greater flexibility 
with dust modelling. Once a data set without extinction 
has been created, multiple samples with different redden- 
ing normalization and modelling can be easily generated. 
Finally, with Galaxia multiple independent random sam- 
ples can be generated, which is crucial for doing a proper 
statistical analysis. Due to the above mentioned advan- 
tages we chose Galaxia as our preferred model to create 
mock samples. 

Table [2] lists the basic parameters for each of the two 
models. For the dust modelling, we chose the default 
value for the Besangon model, where the dust is mod- 
elled by an Einasto disk with a normalization of Ay = 
0.7mag/kpc. This is reasonable for the high latitudes 
that we simulate. Assuming a i?v'=3.1, this corresponds 
to a reddening rate of E{B — V)= 0.23 mag/kpc. No ad- 
ditional dust clouds were added. For the Galaxia model, 
we present results with the dust modelled by an exponen- 
tial disk, with the reddening rate in the solar neighbour- 
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Figure 2. (a) The mock Besangon sample and (b) a mock Galaxia sample for —70 < b < —50, J > 10.3. As in Figure 1, the Aquarius 
stream region is delimited by the red box, with both mock samples displaying a apaucity of stars in this region. A reddening rate of 
E{B — V)=0.23 mag/kpc is used for both the model samples displayed (see Sectionlslfor details). 



Model Solar position Solar motion Vc E{B — V) rate 

(x, y, z) kpc (U, V, VK)kms^"'^ kms^"*^ mag/kpc 

Besangon (-8.5, 0.0, 0.015) (10.30, 6.30, 5.90) 226.40 023 

Galaxia (-8.0, 0.0, 0.015) (11.10, 12.24, 7.25) 226.84 0.23, 0.53, Schlegel 

Table 2 

Parameters for the Galaxia and Besangon models used for comparison with the RAVE sample. 



hood normalized to 0.23 and 0.53 mag/kpc, wh ere the 
latter is taken from Binney & Merrifield (1998). Also, 
we present results for a model where the reddening at 
infinity is matched to that of the value in Schlegel maps. 
To convert E{B — V) to extinction in different photo- 
metric bands we used the conversion factors in Table 6 



of Schlegel et al. (1998). 



3.2. Mock sample generation 

The mock samples were created from Galaxia and Be- 
sangon using analogous methodology. Firstly, to create 
the Besangon sample we queried Al x A6 — 50° x 20° 
regions using the online query form imposing the /-band 
magnitude limits of RAVE of 9 < / < 13, making no 
biases in spectral type. A distance limit oi d = 20 kpc 
is imposed as most RAVE s tars (with the excep tion of a 
few notable LMC stars - see Munari et al. (2009)) should 
be within 15 kpc (Brcddcls ct al. 20TU). (Jrid-steps of 
10° in I and 5° in b were lised m tbe query. 

To generate samples from Galaxia we simply generated 
a full catalog over the area specified by < / < 360, 
—90 < 6 < and 9 < J < 13 and then extracted the 
required samples from it after correcting for extinction. 
Since Galaxia allows oversampling, the initial catalog was 
generated with an oversampling factor of 10, so that later 
on 10 independent random realizations could be created. 

Using Monte Carlo techniques each model was then 
resampled first to create a uniform distribution in I- 
magnitude and then resampled again to exactly mimic 
the shape of the DENIS /-band distribution in a x 
Ab = 50° X 20° region. This ensures that the distance 
distribution will be similar to the RAVE sample. Each 
generated sample is then further reduced to those stars 
with J > 10.3 to mimic our sample selection in Section 
[2) Finally, the number of stars in the mock sample is 



normalized to that of the RAVE sample in sub-regions 
of Al X Ab = 25° X 10°, where this division into sub- 
regions was required to better suit the curved bound- 
ary of the RAVE survey area. For the Besangon sam- 
ple, the / and b co-ordinates were smeared out to remove 
the discretization by adding a uniform randomization of 
the same extent (since the Galaxia sample was already 
smoothly distributed no such procedure was required). 
Also, for Galaxia ten mock data samples were created 
for each dust modelling scenario, enabling a better han- 
dle on the statistical significance of the Aquarius stream. 
Finally, to simulate the RAVE radial velocity measure- 
ment errors a scatter of ct = 2kms^^ was added to the 
models' radial velocities. 

3.3. Statistical significance of Aquarius 

Figure [2] shows the Besangon and one of the Galaxia 
samples (with E{B ~ V) — 0.23 mag/kpc) for the same 
area of the sky as in Figure [l]i. We see that both models 
do a fair job of reproducing the gross features of the data. 
A detailed analysis of comparing both the Besangon and 
Galaxia models to RAVE will be presented in an upcom- 
ing papers by A. Ritter. In this analysis it is sufficient to 
note firstly that the Galaxia model produces a better rep- 
resentation of the density of halo stars (i.e., those stars 
with larger Mos) than the Besangon model. Moreover, 
the Galaxia model better reproduces the Vios distribu- 
tion as a function of Galactic latitudes than Besangon: 
for bins of 25° in Galactic latitude, on average Galaxia 
agrees with the data to within 2kms~"'^ and 3kms~^ for 
mean and dispersion in V|os respectively, compared to 
3kms~^ and 4kms~^ for Besangon. 

To compare the generated samples to the RAVE sam- 
ple, we establish cells of size Al x Al^ios and for each 
cell compare the number of stars from RAVE and the 
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Figure 3. As in Figure [T] but for the latitude 
model with Schlegel et al. ( 1998^ dust mapping, 

the RAVE data m the top panel and are outhned by the green box. The placement of the Aquarius stream from Figure [l] is outlined by 
the red box. Other than the LMC structures are not easily discernible. 



50 < fe < -30 (a, b: top) and -90 < b < 
Schlegel et al. ( 1998^ dust mapping for comparison. LMC stars can be seen at 270° < 



70 (c, d: bottom) using the Galaxia 
< 290°, 230 < Vios < 310kms-i 



mock samples. For each sample in the i-th cell there are 
jyModei g^j^j-g g^jjj estimate the standard deviation by 

fTi = y^]yModoT^ consider an overdensity significant if 



RAVE 



N. 



Model 



(1) 

where A^J^^ve ^^.^ number of RAVE stars in the i-th 
cell and iVf °doi jg ^^^^ler N^"^ or ivf where g = 1..10 
in the latter signifies the sample number from Gala xia. 
Following a procedure similar to Helmi et al. (19991, we 



identify overdense regions in the Galactic latitude slice 
—50 < b < —70 by varying the cell sizes with longi- 
tude slices ranging from Al = 25, 35, 50, 70 and ra- 
dial velocity bins ranging from AVios = 20, 25, 35. ..100. 
We then evaluate the percentage of the various cell 
sizes which identify the region around 30° < I < 75°, 
—250 < Vios < — 150kms~^ as having a 4(t deviation. 
As we have 10 samples for Galaxia, we take the aver- 
age over all the samples, obtaining a mean and standard 
deviation for this value. 

The following results are found: using the Besangon 
model, 96% of the different cell sizes identify that the 
number of stars in the data are Aa overdense around 
Aquarius compared to the model. For Galaxia using 
E{B -V) = 0.23 mag/kpc, we find that 80 ± 15% of ceU 
sizes give Aquarius as a 4cr deviation, while E{B — V) = 



0.53 mag/kpc gives 75 ± 15% and the Schlegel results 
yield 80 ± 18%. How the dust is modelled at these high 
Galactic latitudes therefore has little impact on the re- 
sults. In general we can conclude that the models ro- 
bustly show that there is a statistically significant con- 
centration of stars at the Aquarius stream's location. In- 
deed, for some cell sizes and models the overdensity can 
be as high as 11a. This confirms what can be seen by 
eye: the stream as an overdensity in the outlying regions 
of the velocity distribution. 

3.4. Localization of the stream 

In addition to identifying the statistical significance of 
the Aquarius stream, we also used the models to search 
for additional members of the stream and possible re- 
lated substructures. We compared the RAVE data and 
mock Galaxia samples for surrounding latitude cuts of 
-50° < 5 < -30° and -90° < 6 < -70°, where once 
again we consider only those stars with J > 10.3. Figure 
[3] displays the data in these two latitude ranges com- 
pared to Galaxia models using the Sch legel dust model. 
We repeated the analysis of Section [373] looking for over- 
densities in the RAVE data compared to the Galaxia 
models, varying the cell size and identifying regions with 
repeated 4cr signals. 

For the —50° < b < —30° sample, the region around 
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Figure 4. Le/t; MDF for the Aquarius stream members (red line), whose typical metallicity uncertainty is A([M/H]) ~ 0.2 dex (ztlcr 
shown). The MDF of other halo stars with -70° < < -50°, J > 10.3, |Vios| > 200 kms'^ is shown for comparison (dotted line). Right: 
Teff - logg plane for RAVE stars in the region —70° < fe < —50°, 30° < I < 75°, J > 10.3. Stream candidates are highlighted as solid 
red points and a Padova isochronc with lOGyr, [M/H] = —1 overplotted. The yellow region indicates la in both T^g and logg from this 
isochrone. 



270° < / < 290°, 230 < Mos < 310 km s"^ is consistently 
identified for all the various dust models as significantly 
overdense: on average 95% of cell sizes identify this re- 
gion as containing a 4cr signal. These stars are associated 
with the Large Magellanic Cloud (LMC) and it is reas- 
suring that our technique detects it. 

For both latitude ranges, -50° < b < -30° and 
—90° < b < —70°, there are no detections of statisti- 
cally significant overdensities in the vicinity of the Aquar- 
ius stream's velocity and longitude range. Also, for the 
—90° < b < —70° sample no particular region has consis- 
tent 4(T deviations when compared to the Galaxia models. 
The region around -50° < 6 < -30°, 320° < / < 350°, 
150 < Vios < 300kms~^ is detected in - 50% of the 
trials as being overdense for this latitude cut, irrespec- 
tive of dust modelling. A similar detection is also found 
for the region -70 < b < -50, 260° < I < 340°, 
100 < Vios < 300kms"\ in the same latitude range as 
the Aquarius stream. These detections are not as sig- 
nificant as Aquarius and are in a different region of the 
sky. 

In general, we find that there are no stars clearly associ- 
ated to the Aquarius stream in adjacent latitude cuts; no 
overdensities were detected in the vicinity of the stream's 
velocity and longitude. This may be caused in part by 
the survey boundary, but the sharp localization of the 
stream is nevertheless intriguing. In Section [6] we fur- 
ther investigate the localization of the Aquarius stream, 
examining its possible relation to the two marginal over- 
densities detected above. 

4. POPULATION PROPERTIES OF THE 
AQUARIUS STREAM 

4.1. Metallicity and logg-T^ff plane 

RAVE gives estimates of stellar parameters from the 
spectra which we can use to establish the basic proper- 
ties of the population of the Aquarius stream. Conserva- 
tively the stellar parameters are accurate to ~ 0.2 dex in 
[M/H], 400 K in Tes and 0.5 dex in logg when compared 
to external measurement s, though internally the errors 
are significantly smaller (Zwitter et al. 2008). For 13 



of the 16 stream candidates we have estimates of stellar 
parameters. 

Figure [4] shows the metallicity distribution function 
(MDF) and the log g-Tes plane of these stars where we 
compare the latter to the background population. The 
conservative estimates of the errors are also shown. Note 
that we do not apply the metallicity calibration of Equa- 
tion 20 in [Zwitter et al. ( 2008 1 as this calibration does 
not extend down to halo metallicities [M/H] < —1.5. 
Therefore, the derived MDF is best seen relative to 
background halo stars. We plot the MDF for stars 
selected -70° < b < -50°, ^ < 30° or Z > 75°, 
J > 10.3, \Vios\ > 200kms"^ The stream's MDF 
peaks at a slightly higher metallicity than these back- 
ground halo stars with a slightly tighter distribution: the 
stream has an average [M/H] = — 1 ± 0.4 compared to 
[M/H] = -1.1 ± 0.6 for the background. Both distribu- 
tions show metallicities with are rather high for the halo. 
The RAVE 3rd data release, Siebert (2010, in prepa- 
ration) show that these RAVE stellar parameters tend 
to overestimate the metallicity of stars with low SNR, 
and effect that would be on the order of 0.1 dex for the 
stream stars. This data release will present improved 
stellar parameter from a modified pipeline, as well as 
a new metallicity calibration from an extended metal- 
licity range. Hence, these results should a better han- 
dle on the stream's MDF. Clearly, however, follow-up 
high-resolution spectroscopy is required to derive accu- 
rate abundances to better understand the group's chem- 
ical abundance properties. Nevertheless, from the initial 
RAVE metallicities, we can conclude that the stream's 
MDF is consistent with bac kground halo stars. 

Using Padova isochrones ( Girardi et al.|[2002 ) we find 
the best fitting isochrone to be that of lO Cyr, [M/H] = 
— 1 which is over-plotted in Figure |4] (right), as well 
as a highlighted region showing the ±la bounds in 
(log 5, Tcff) from this curve. Most of the stream stars 
fall within this region, though clearly the isochrone fit 
is preliminary given the size of the stellar parameter er- 
rors. From both the MDF and the isochrone fit, however, 
there is a general indication that the Aquarius stream is 
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metal-poor and old. 

4.2. Isochrone derived distances 

We use the isochrone fit from above to derive distances 
to the candidate stars, using the J-band magnitude. To 
derive Mj from the isochrone we find the nearest point 
along the isochrone to the actual data point by mini- 
mizing the distance in log g and Tcff between them, nor- 
malized by the standard error in each. Extinction is of 
the order of Aj = 0.04 ± 0.01, and is calculated iter- 



atively from the distances using Schlegel et al. (19981 
d ust maps and assuming a Galactic dust distribution as 



Beers et al. 



(20001. Errors are calculated via Monte 
Carlo, generating lUO points around the data point with 
(TTeft = 400 K and ciogg = 0.5 dex, propagating through 
to a distribution of distances, from which a standard de- 
viation is derived. 

The distances are listed in Table [3] as dj, where the 
distances range from 0.4 to 9.4 kpc (distance moduli: 
m — M ~ 8.3 to TO — M ~ 14.9), with a mean distance 
of dav = 3.8 ± 3.2 kpc {m ~ M ^ 12.9). There is a hint 
of a bimodal population of closer (sub- and red clump 
giants) and farther stars (tip of the giant branch). How- 
ever, given the uncertainties in the stellar parameters 
these distances are uncertain and the reality of this bi- 
modality is therefore debatable; in the next section we 
develop another distance estimate which has a smoother 
distribution function. 

The large distance range raises the question whether 
the Aquarius stream is a distinct entity or comprised 
of multiple structures. The high-resolution abundances 
mentioned above would help answer this question by as- 
certaining if the group has any distinctive chemical sig- 
natures compared to other halo stars. Occam's razor 
would weigh against two structures forming this local- 
ized stream however. Further, in Section [6.1| we develop 
a model for the Aquarius stream under the assumption 
of a single satellite dissolving in the Galaxy's potential. 
The model predicts that the stream is spread in XYZ 
away from the sun, with distance dmodci = 3.2 ± 0.8 kpc 
in the direction 30° < / < 75°, -70° < 6 < -50°. The 
distance range derived above therefore probably reflects 
more on the distance errors than the real distribution for 
the stream. We assume that the Aquarius stream is a 
single, distinct object. 

The isochrone from Figure |4] has a /-band turn-off of 
Mj = 3.5. Hence, for the distance moduli above we 
could expect turn-off stars in the range / = 11.8 - 18.5. 
The lower magnitude falls within the RAVE magnitude 
limits (9 < / < 13). However, RAVE's unbiased selection 
criteria mean that the thin disk dominates dwarf /turn- 
off stars, even at these higher magnitudes. Our sample 
of halo dwarfs is therefore too small to detect the turn- 
off, and we only see giant stars in our Aquarius stream 
sample from RAVE. 

4.3. Reduced Proper Motion Diagram 

The coherence of the group selection is shown by the re- 
duced proper motion diagram (RPMD), which plots the 
reduced pr oper motion (RPM) a gainst color. Described 
in detail in lSeabroke et al.| ( [2008l ), the RPMD essentially 
creates a MR diagram trorn the proper motions, where 
the absolute magnitude is smeared by the variation in 




Figure 5. Reduced proper motion diagram for tiie background 
RAVE stars (blaclc points) and the Aquarius stream stars (red 
points). Tlie isoclirone from Figure |4] is plotted with a tangential 
velocity of vt = 250kms~^ (solid line), vt = 150kms~^ (dotted 
line) and = 350kms~^ (dashed line). The coherency of the 
group is clear also in this diagram. 



the tangential speed of the stars. Halo stars have a large 
dispersion in tangental velocity and so this smearing is 
large. In contrast, for a small, nearby section of a stream 
the transverse velocity spread is small and we effectively 
recover magnitudes for the stars. The RPM is given by 

i/j = J + 5 log ^ + 5 = Mj + 5 log Vt - 3.379 (2) 

where J and Mj are the apparent and absolute mag- 
nitudes respectively, /x is the proper motion in arcsec 
yr~^, Vt is the tangential velocity in kms""'^. Here we 
have again used 2MASS colors. Thus, from the observ- 
ables J and /i we can establish something about the more 
fundamental parameters Mj and vt without requiring a 
either distance or a radial velocity. Figure [5] gives the 
RPMD for the stars in our magnitude and latitude se- 
lected sample with the Aquarius stream candidates over- 
plotted, where for the latter the more accurate PPMXL 
proper motions were used. Note that for the distances 
of these stars the reddening will also be of the order of 
E{J — K) ^ 0.02, which does not effect the plot sig- 
nificantly and is neglected. The isochrone from Figure 
|4] is over-plotted, where we find that a large tangential 
velocity of vt — 150kms~^ to 350kms~^ is required to 
shift the isochrone to a reasonable fit, which compare to 
Vt 230±100kms"^ for the halo for (l,b)=(55°, -60°). 
Once again, the group is consistent with a halo stream. 

We will see in Section [6] that the tangential velocity 
for the stream is indeed within a relatively narrow range 
as shown in Figure [5j A few of the redder, more-distant 
giants deviate from the rest of the group but they also 
have larger errors in their proper motions, which trans- 
lates into larger RPM errors as shown: they are within 
2cr of the group fit. The consistency of the fit for the 
bluer (nearer) stars supports their inclusion in the candi- 
date list, though the range in values for the RPM of halo 
stars means that we cannot exclude contamination from 
the halo. Indeed, from Figure[l]it is clear that we expect 
a few of the stars to be non-members. We therefore take 
the consistency of the RPMD to be a good indication 
of the consistency of the Aquarius member selection but 
not absolute proof of membership. 
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ID 


Obadatc 


J 


K 


T^ff 
K 


[M/H] 


log 9 


dj 

kpc 


dp{ 
kpc 


kpc 


d2 
kpc 




^BB 
kpc 




J221821. 2-183424 


20060602 


10.34 


9.68 


4572 


-1 


54 


1.06 


5.6 ±2.1 


2.4 ± 1.0 


9.6 ± 2-1 


6.7 ± 2 


9 


2.8 ± 


5 


C2222241-094912 


20030617 


10.64 


9.79 












16.3 ± 11.9 












C2225316-145437 


20040628 


10.34 


9.57 


4104 


-1 


29 


1.01 


7.3 ± 2.7 


3.4 ± 1.5 


10.3 ± 2.4 


7.9 ± 3 


1 


1.4 ± 


1 


C2233207-090021 


20030617 


11.66 


11.28 












0.8 ± 0.3 












C2234420-082649 


20050914 


10.67 


10.13 


5263 


_2 


02 


2.43 


1.9 ± 0.7 


1.9 ± 0.8 




3.1 ± 2 


8 






J223504. 3-152834 


20060624 


10.36 


9.65 


4795 


-0 


33 


3.05 


1.0 ± 0.5 


2.5 ± 1.1 




1.3 ± 1 









C2238028-051612 


20050807 


11.53 


10.74 


4606 


-0 


86 


1.49 


7.1 ± 2.7 


6.7 ± 5.4 












J223811. 4-104126 


20060804 


10.42 


9.90 


5502 


-0 


78 


4.16 


0.4 ±0.1 


0.5 ± 0.2 


0.3 ± 0.1 


0.5 ± 


9 


0.4 ± 


3 


C2242408-024953 


20050909 


11.63 


10.82 


4159 


-0 


75 


1.53 


9.4 ±3.4 


13.4 ± 8.7 




8.5 ± 5 


6 






C2246264-043107 


20050807 


11.26 


10.72 


5142 


-1 


22 


2.65 


1.8 ± 0.9 


1.8 ± 0.7 


3.9 ± 2-1 


3.2 ± 3 





2.0 ± 


8 


C2306265-085103 


20030907 


10.31 


9.47 












3.3 ± 1.5 












C2309161-120812 


20040627 


10.68 


9.97 


5219 


-0 


66 


2.94 


1.0 ± 0.6 


5.5 ± 2.9 




1.5 ± 1 


4 






C2322499-135351 


20040627 


10.82 


10.28 


5043 


-0 


64 


2.45 


1.9 ± 0.9 


2.1 ± 0.8 




2.3 ± 1 


6 






J232320. 8-080925 


20060915 


10.96 


10.47 


5286 


-1 


10 


3.50 


0.7 ± 0.4 


0.8 ± 0.3 


1.0 ± 0.5 


1.8 ± 2 


2 


1.1 ± 


7 


J232619. 4-080808 


20060915 


10.51 


9.76 


4225 


-1 


22 


1.14 


6.7 ± 2.7 


5.6 ± 2.9 


10.6 ± 2.8 


7.1 ± 3 





3.5 ± 


6 



Table 3 

The Aquarius Stream candidates selecte d fr om the RAVE data and their parameters. The dis tanc es dj are derived fro m the isoc hrone in 
Figure B| while cLr are derived in Section | 4.3[ The extra distances are dn lIBr eddels et al.|2010[l. di llZwitter et al.|2010|l and JBurnettl 
I fc Binney|20'rny : 



4.4. RPM derived distances 

If we accept the group's tangential velocity of vt = 
250kms~^, we can use the RPM to establish a second 
estimate of the distance to the stars. From Equation [2] 
we have the distance modulus 



J ~ Mj ^ J - Hj + 5\ogVT - 3.379 



(3) 



From this we have distance moduli ranging from 8.5 to 16 
for the group members. The corresponding distances are 
listed in Table [3] as dn , with the errors calculated using 
the upper and lower tangential velocity bounds as well as 
the proper motion errors in Hj. These distances differ 
somew hat f rom those calculated using the isochrones in 
Section [42] but are of the same order of magnitude, with 
an average value of dav = 4.5 ± 4.6 kpc (/i — 13.2). Two 
stars (C2222241-094912 and C2242408-024953) have very 
large distances but the errors are also large. If we exclude 
the largest of these but retain the other for consistency 
with the isochrone distance average, the mean distance 
reduces to dav = 3.6 ± 3.4 kpc (/i = 12.8), which is very 
similar to the value found for the isochrones. 

4.5. Comparison with other distances 
In Table m we also list distances derived in[Brcddcls et 



al] IgOl Ol (dB) jZwitter et al.| ( |2010| (dz) and |Burnett 
& Binney (2010]) (dBs), where tnese distances are all de 



rived from RAVE stellar parameters, employing various 
methodology. Comp aring the above di and dpt to the dis- 



tances calculated in Zwitter et al. (20101, for which 11 
of the 15 Aquarius stars have an entry, we find that the 
isochrone distances agree better with a ± cr for the dif- 
ference (dz-di)/dz of 23% ±20%, while for {dz-dn)/dz 
we have —10% ± 100%. This is somewhat unsurpris- 
ing given that the Zwitter distances and the isochrone 
distances are both based on RAVE stellar parameters. 
Interestingly, however, for 6 stars that have distance s 
calculated by the method of Burnett & Binney (2010), 
the RPM distances fare betteE (deB — di)/dBB gives 
-100%±170% whi le (dBB- dR)/dRR yields -30%±70% 



For the 6 Breddels et al. ( [2010) entires we have much 
larger discrepancies ot (dB--di)/dB of —200% ±40% and 
(ds — dR)/dB of 330% ± 60%. Clearly, all these discrep- 
ancies imply that the individual distances listed in Ta- 
ble [3] have large uncertainties. In general, however, the 
RPM distances give more consistent kinematics than the 



isochrone distances as we will see below in Section [S] 

5. POSSIBLY RELATED SUBSTRUCTURE 

In this section we seek connections between the Aquar- 
ius stream and other known kinematic and spatial sub- 
structures nearby in the Galaxy. We start with the spa- 
tially detected substructures before returning to kine- 
matically detected solar neighbourhood features. 

5.1. Large stellar streams and features 

The nearest companion of the Milky Way, the Sgr dSph 
(Ibata, Gilmore & Irwin 1994), has shed significant de- 
bris on its polar orbit aro und our Galaxy. The all-sky 
mapping of this debris by Majewski et al. (2003) using 
2MASS M-giants clearly showed the pla ne or the Sgr de- 
bris. F urther studies such as t hose by [Newberg et al.| 



have revealed further branches and details withm the de- 



(| 2003| ; [Martmez et al. ( |2004[ ); [Belokurov et al.| p"006 l 

bhe 



bris wraps. Recently, Yanny et al. (2009a) used M and 
K gia, n ts selec ted from SDSS and SEGUE data |Yanny| 
et al. (2009b) to provide additional observational con- 
straints on the stream. 

The Aquarius stars fall fairly close to the orbital plane 
of the Sgr dwarf. Also, the isochrone fit from Section |4] 
is consistent with a p opulation of lOGyr, [M/H] = — 1. 
Layden fc Sarajedini| |2000| ) obtain CMDs of Sgr field 
populations, finding a dominant old and intermediate age 
popul ation of llGyr, JM/H] = -1.3 and 5 Gyr, [M/H] = 
—0.7. Giuffrida et al.| (2010l find a range of populations 
in the periphery ot Sagittarius, with [M/H] = —2.34 to 
-0.6, while the dominant population has a similar metal- 
ficity to 47 Tuc with [M/H] = -0.6. Given the errors, 
the isochrone fit for the Aquarius stream is consistent 
with the Sgr dwarf. We thus investigated a possible link 
between the Aquarius stream and the Sagittarius dwarf 
debris. The details of this investigation are given in Ap- 
pendix A. 

The overall result is that the Aquarius stream's kine- 
matics do not match those of the Sagittarius dwarf 
debris, calculated using a variety of potential models 
(oblate, spheroid, prolate, triaxial) from Law (2005, 
2009). The oblate model shows a potential match for a 
small section of nearby debris when considering the line- 
of-sight velocity in the Galactic rest-frame, V^ai, alone. 
However, the full kinematics of V^, Vr, Vz displays that 
the kinematics of the Aquarius stream and this nearby 
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Figure 6. (a) L^-Lpcrp and (b) L^-Energy (Lindblad) plane for the Aquarius stars with the solid red points calculated using and the 
open red points dj. In the background the mostly thin-disk GCS stars are shown as black points. The uncertainty in the solid red points 
are shown via the clouds of small colored dots, which give the Icr spread for the MC simulation. Each color shows the spread for a single 
solid red point. A high degree of co-variance in the large errors is evident, though the retrograde nature of the Aquarius stream stars is 
clear. The solid curves in the Lindblad diagram represent the circular orbit loci for this potential. 



section are actually quite differenlp^ The possible con- 
nection is further ruled out by theiact that the oblate 
halo potential model does not compare well with other 
observational data for Sgr dwarf debris. 

Since the Aquarius stream lies in the southern part 
of the RAVE data, it could not be discovered in the 
main, northern SDSS survey. Thus the stream is far re- 
moved from the SDSS-discovered substructures, includ- 
ing the Canis Major ov erdensity at {I, b) = (240°, -8°) 

and the Virgo overdensity at 



( Martinez et al. 
(l, b) = (300°, 



2005|) 



2008[ ). Further, the 



(Juric et al. 

stream is located between the southern' SEGUE SDSS 
stripes so it unsurprising that this has not been de- 
tected in this survey. The stream's Galactic latitude of 
b = —60° rules out a rel ation to the more plana r Mono- 



ceros stream {b < 40°) ( Penarrubia et al.]|2005 ). Its ve- 
locities and latitude are also inconsistent with the thick 



disk asymmetries detected by Parker et al. ( 2003 , 2004 ) 



The Hercules- Aquila cloud, again detected usmg SUSS 
photometry, is located at I — 40 



below the plane by 50° Belokurov et al 
velocities of the b > 



and extends above and 

The 
1 



(2007) 



0° segment are V^ai — +180 km s 
and the structure ranges over heliocentric distances of 
d — 10 — 20 kpc. The Hercules- Aquila cloud is near the 
Aquarius stream on the sky. However, despite the lack 
of velocity data below the plane, it can be clearly seen 
that the two entities are separate: the centering in (1, b) 
for the two are shifted from each other and their distance 
ran ges are clearly incompatible. Additionally, in Section 
|6.1| we trace the orbit of a simple model for the Aquarius 
stream and the resulting region of phase-space that the 
debris inhabits does not overlap with the Hercules- Aquila 
cloud in {I, b, Vgai). 

5.2. Solar neighbourhood streams 

We have calculated orbits for candidate stars in the 
potential of Helmi et al (2006), which has contributions 
from a disk, bulge, and dark halo. Table |4] gives averages 

We use the [Pehnen &: Binney| ( |l998| l values for the solar pe- 
culiar velocity of {U, V, W) = 10, 5, 7kms~^ with respect to the 
LSR, which we set at a rotation velocity of 220kms~^. 



for various quantities derived from these orbits as well as 
the median quantities for the overall kinematics, using 
both sets of distances. Note that we chose the median 
as it gives more consistent results, and for this reason 
we also excluded the two most distant stars with d > 9 
kpc as their kinematics differed greatly from the others. 
Also, the values for the pericentre and apocentre only 
include non-radial orbits. 

Figure [6] shows the L^-Lpcrp and L^-Energy (Lindblad) 
planes for orbits based on both distance estimates, where 
to aid compariso n to other stud i es we use here energies 
as calculated in Dinescu et al. (1999). Note that the 



scatter of the isochrone distance results is large so the 
majority of these points lie off the plot, as do some of 
the RPM distance results. We p lot for reference stars in 
the Geneva Copenhagen survey ( Nordstrom et al.|2004 ), 
which is comprised mainly of thin and some thick disk 
stars. The circular orbit loci for this potential are also 
shown in the Lindblad diagram. To show the typical 
error covariance we also ran a Monte Carlo (MC) simu- 
lation for each star (with the RPM distances). From the 
errors in distances, proper motion and radial velocity we 
generated a sample of 1000 points representative of each 
distribution, which were then propagat ed through into 
the momenta and energy. Following Wylie-de Boer et] 
al. (2010), in Figure [6] we plot the resulting distributions 



within 1(7 of each ofthe average values. These demon- 
strate the large non-Gaussian errors and covariances in 
Lz, Lpcrp, Energy: the Aquarius stream could not be 
found initially in these planes. Indeed, all three values 
are quite ill-constrained with the current uncertainties in 
the stellar distances. 

Nevertheless, we can at least say that the stars are 
retrograde and that they ar e away from the notable 
halo feature of Helmi et al. (1999) at {L^, ipcrp) = 

rj 



(1000, 2000) kpc km s ^. The Aquarius stream is also 
not near th e prograde and retrograde features of Kepley 



et al. 



(2007) at (L,, Lp 



(2500, 500)kpckms"' and 



{Lz, L 



porp; 



(-2500, 1700) kpc km s"'. Also, with val- 
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Table 4 

Orbital properties of the Aquarius stream using (A) isochrone distances, (B) RPM distances and (C) the satellite model in Section |6.1 

selected between -70 < fe < -50, d < 5 kpc, ts = -700 Myr. 



ues of (Vaz, Vae, v) = (125kms-\ 205kms-\ -15 _ 
it is not one of the newly d etected solar ne ighbourhoo 



streams listed in Table 2 of Klement et al. (2009). An- 



other solar neighbourhood stream i s the Kapteyn group 



which Wylie-de Boer et al. 


(2010D su 


Kgests is stripped 


from uj CJentauri (also see Meza et al. 


(20051), which in 



dwarf galaxy ( Bekki fc Freeman|2003 ). Wylie-de Boer et 



al. ( 2010, ) employ the same potential as we do here, tind 



ingi^ = -413kpckms"\ £; = -1.3xl0'^ 



km^ s ^ for uj 

CenandL^ 200kpckms"\ E 9.5 x 10^ km^ s^^ 

for the Kapteyn group. The Aquarius stream is some- 
what similar to the distribution in L^-Energy for the 
Kapteyn group/w Cen. However, in Section 6.2 we will 



see that our model for the stream rules against an as- 
sociation. Another sig nificant halo clumping found by 
Majewski et al. (19961 towards the north Galactic pole 

— 55kms~^, which is 
The mean Izl for this 



has a retrograde orbit with ^ 
consistent with that of Aquarius, 
moving group of |zav| — 4.5 kpc is rather high when com- 
pared to (|zniedian|, |^av|) ~ (2, 3.5) kpc for Aquarius, as 
well as the Zmax values in Table |4j Moreo ver, our model 
for the Aquarius stream in Section |6.2| does not over- 
lap with the north Galactic pole, again ruling against an 
association. 

6. NATURE OF THE STREAM 

The distances calculated for Aquarius stream stars 
place it fairly close to the sun. If they are of the cor- 
rect order of magnitude then we could possibly expect 
additional stream members in other areas of the sky. 
However, our expl orat ion of the two bounding latitude 
ranges in Section |3.4| yielded no striking overdensities 
that we can immediately associate with the Aquarius 
stream. Two other regions had overdensities detected 
for J > 10.3, though they are not as conspicuous as 
Aquarius: the region around -50° < 6 < -30°, 330° < 
/ < 345°, 190 < Vios < 270 km s"^ (Region A) and 
the region -70° < 5 < -50°, 280° < I < 315°, 
40 < V|os < 130kms~^ (Region B). To establish if these 
additional areas could be associated with the Aquarius 
stream, and if not, how the stream's localization arises, 
we created a simple model of a satellite dissolving in the 
potential of the Galaxy. 

6.1. Model satellite disruption 



Klement et al. 



1 2009 I define V, 

Vae = V^^^ + 2(y VlsrP and u = 



z = V(V + VtsR)2 + H/2, 
arctan((y-f yLSR)/W^) 



To generate a simple satellite dissolution, we first chose 
one of the average, stable orbits - that using d^. for the 
star C2322499-135351 - and integrated the orbit back in 
time. Centering on the orbital positions at various times 
in the past, tg, we generated lO'' test particles from a 
Gaussian sphere with core radius and internal velocity 
dispersion Vc — 300 pc, ay = lOkms"^. Neglecting 
self-gravity we then integrated the orbits of the satellite 
forward in time until the present day. This approximate 
approach suffices as our aim here is illustrative rather 
than finding the definitive orbit for the Aquarius stream. 

Figure [7] shows the distribution of particles for two 
different starting times: one starting two orbital pe- 
riods ago, ts = —700 Myr, and another starting at 
ts = — 5Gyr. These were selected to illustrate two dif- 
ferent extremes, one in which the stream at the present 
day has yet to be significantly phase mixed and the other 
when it is completely phase mixed. In both scenarios the 
Aquarius stream stars sample a volume near the apoc- 
enter of the orbits. The furthermost stars have Z values 
that are substantially larger than Zmax. = 4 kpc for the 
dissolving satellite, a discrepancy that may be resolved 
either by a more detailed model or by more accurate dis- 
tances. In general, however the majority of Aquarius 
(and RAVE) stars are within a few kpc of the sun, and 
only a small portion of the total volume traced by the 
orbits falls within this sample volume. 

In Figure [8] we plot l-b and l-Vios planes for test par- 
ticles within d < 5 kpc of the sun for the simulations 
with the two different starting times, where we chose this 
distance limit as it encompasses ~ 70% of the Aquarius 
stars. We also plot the Aquarius candidates and the stars 
that fall within Regions A and B, as well as the back- 
ground population of RAVE stars. The distribution of 
particles for the simulation with tg = —700 Myr occupy 
a small region in l-b-Vios that does not coincide with the 
locations of Regions A and B. Furthermore, the distri- 
bution in /-Vios remarkably mimics that of the Aquarius 
stream and a picture emerges of how the stream can be so 
localized: with the RAVE survey data we miss portions 
of the stream due to the location of the survey boundary 
in some regions and in others, because Vios overlaps with 
that of the main distribution so the stream is difficult to 
detect. 

In the phase-mixed {ts — ~5 Gyr) scenario we see that 
there is overlap between the region in l-b plane occupied 
by the test particles and for Regions A and B, though 
the V\os values for Region B agree better with those of 
the test particles than Region A. However, it would be 
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Figure 7. XY (top) and XZ (bottom) planes of a simple satellite model of the Aquarius stream at the present day with two different 
starting times: tg = —700 Myr (a, c: left) and tg = — 5Gyr (b, d; right). The position of the sun is marked as a yellow point and the 
Aquarius stars (using d^) as red points. 



difficult to associate both Regions A and B with Aquarius 
without others regions being populated. In particular, we 
could expect overdensities at -50° <b < -30°, / > 240°, 
Mos = +200knis~^ as well as a population at —30° < 
b < -50°, I < 60° out to Vios = -200knis^\ Such 
overdensities are not observed in the data. 

These simple models of a dissolving satellite therefore 
suggest that the localization of Aquarius is due to further 
regions of phase-space not yet being populated: the re- 
gion in {I, b, Vios) space occupied by the tg = —700 Myr 
simulation is more consistent with the observed popula- 
tion than that of the tg = —5 Gyr simulation. This sug- 
gests that the stream is most likely dynamically young, 
resulting from a recent disruption of a progenitor, and 
has not yet undergone phase mixing. Indeed, our simple 
model of a recently disrupted satellite is very successful 
in reproducing the main features of the Aquarius stream. 
In Table |4] we therefore list the parameters found for test 
particles from this model in the Aquarius stream's lati- 
tude range. These are quite similar to those found using 
the di and dji. The results also co rrob orate our obser- 
vation from the RPMD in Section 14.31 that the stream 
stars have a const ant t ransverse velocity, a fact which was 
used in in Section [T4] to derive d^ . Test particles for the 
tg — —700 Myr simulation within the Aquarius stream 
latitude range and distance range, i.e., —70° < b < —50° 



and d < 5 kpc have Vt = 250 ± 33 km s ^ . This agrees 
with the value of vt = 250 ± 100 found from the RPMD. 

From Figure |8] we see that the recent-disruption model 
suggests that at -10° < 5 < -^30°, 330° < I we can 
expect a smaller population of stars associated with the 
Aquarius stream out to Vios = -t-350kms~^. This begins 
to overlap at -1-20° < b with the RAVE survey area, 
though we do not detect such a population in the data. 
Future releases of RAVE data with more observations in 
this area, combined with more careful modelling of the 
stream, will enable a better understanding of whether 
this area is indeed populated by Aquarius stream stars. 

6.2. Progenitor of Aquarius 

The above scenario of a dynamically young stream is 
not inconsistent with an age of 10 Gyr for the Aquarius 
candidates, as estimated from the isochrone fit in Section 
14.11 the stream can be seen as a remnant of an old satel- 
lite that has been recently disrupted. As to the nature 
of the progenitor of the Aquarius stream, it could either 
be a dwarf galaxy or a globular cluster. The survival 
of this progenitor would depend on its concentration: it 
could either have been tidally stripped or have undergone 
complete disruption. 

To search for possible globular clusters that the Aquar- 
ius stream could have been tidally stripped from, we per- 

( |l996| ) catalog of known 



formed a search of the Harris 
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Figure 8. The l-b (top) and l-Vioa {bottom) planes for the present-day satellite particles within d < 5 kpc for simulations with two 
starting times, ts = —700 Myr (a, c: left) and ts = —5 Gyr (b, d: right). The test particles are the small colored points with the color 
denoting their Galactic latitude; those at b = +20° are green, progressing to light blue at 6 = —10°, dark blue at fe = —40° and violet 
at 6 = —70°. The Aquarius stream stars (large red points) are superimposed as are the two other regions with possible overdensities: 
-50° < b < -30°, 320° < I < 350°, 150 < Vios < 300kms-i (Region A: black squares) and -70° < b < -50°, 260° < I < 340°, 
100 < Vios < 300kms-l (Region B: small orange diamonds). RAVE stars having J > 10.3 are plotted in the background (grey points). 
The Aquarius stream bears a striking resemblance to the left-hand scenario, suggesting the stream is not yet phase mixed. 



globular clusters, selecting those with —1.5 < [Fe/H] < 
—0.5, 1.5 < i?GC < 9 kpc and Z < 4 kpc, where the 
latter limits are taken from the model satellite orbit in 
Table |4| We then compared the distribution of the clus- 
ters in /, b, V\os to that of the model satellite stream 
{tg — —700 Myr), and found no globular clusters that 
match the simulation. Also, uj Cen, with an I, b, Vios 
of (309°, +15°, 233kms"^) and = -413kpckms~\ 
is not consistent with not-yet-phased-mixed scenario in 
Figure [8] and Table |4] With our uncalibrated metallici- 
ties it is difficult to compare the MDF to that of w Cen. 
A high- resolution spectrosco pic abundance study, such 



Wylie-de Boer et al. (2010), is required (as well as fur- 



ther modelling of Aquarius) to definitively understand if 
Aquarius is related to lo Cen. 

The progenitor of the Aquarius stream therefore is cur- 
rently unknown. 

7. CONCLUSION 

In this paper we report the detection of a new halo 
stream found as an overdensity of stars with large he- 
liocentric radial velocities in the RAVE data-set. The 
detection is enabled by RAVE's selection criteria creat- 



ing no kinematic biases. The fifteen member stars de- 
tected have Vios — —199 ± 27kms~^ and lie between 
-70° < & < -50°, 30° < / < 75°, J > 10.3 in the 
constellation of Aquarius. We established the statisti- 
cal significance of the stream by comparing the RAVE 
data, in the Galactic latitude range —70° < 6 < —50°, 
to equivalent mock samples of stars created using the 
Besangon Galaxy model and the code Galaxia. For dif- 
ferent cell sizes, Al x A Vios, we compared the number of 
stars in the data and models, finding that for the major- 
ity of cell sizes the region around the Aquarius stream 
exhibited a 4tT overdensity in the data, irrespective of the 
dust modelling. Searching for additional overdensities in 
neighboring latitude regions yields no structures of the 
same level of significance (other than the LMC), though 
two regions are identified as being marginally overdense. 

For most of the Aquarius stars RAVE stellar param- 
eter estimates are also available. The member stars are 
metal-poor with [M/H] = — 1 ± 0.4 and we derive a 
preliminary isochrone fit in the Toff -log. 9 plane with an 
population age of 10 Gyr. Both the Vios and metallic- 
ity are consistent with the group being within the stellar 
halo. We further use a Reduced Proper Motion Diagram 
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to derive the transverse velocity for the stream, finding 
vt = 250±100knis^^ for the group. This again places it 
within the Galaxy's halo. We use the isochrone fits and 
the RPMD to provide distance estimates to the stars, 
where we prefer the latter as they give more consistent 
kinematics. 

We investigated the relation of the stream to known 
substructures. We first discussed the probability of the 
stream being with debris from the Sagittarius dwarf. 
This is a priori plausible because the stream does not 
fall far from the orbital plane of the Sgr dwarf and the 
stream's metallicity is consistent with that of the dwarf. 
A co mparison to the mod els of Sagittarius dwarf debris 
from Law et al. (2005) and Law et al. (2009), shows that 
although the majority of models do not yield a good fit, a 
certain selection of nearby stars in the oblate model pro- 
vides a reasonable fit in the AQ-V^gai plane. This is most 
likely just coincidental however: the distributions in both 
distance and are clearly inconsistent with those of the 
Sgr stream. Also, the oblate model is the least favoured 
of all the models when compared to the most recent data 
for the Sagittarius stream. We thus conclude that the 
Aquarius stream is most likely not associated with the 
Sagittarius dwarf. A search of other known substructures 
both in the solar neighbourhood (e.g. Kapteyn group) 
and in the solar suburb (e.g. Canis Major and Virgo 
overdensities, Monoceros stream, Hercules- Aquila cloud) 
yielded no positive identifications. 

Finally, to understand better how the stream is both 
local and localized on the sky, we performed simple dy- 
namical simulations of a model satellite galaxy dissolv- 
ing in the Galactic potential. We presented simulations 
for two time-scales, one where the satellite is dissolving 
and the other when it is completely phase mixed. We 
compared the distribution in 6, Vios space of nearby 
tracer particles at the present day to that of the Aquar- 
ius stream stars plus the two other marginally overdense 
regions found in the RAVE data. The model in which the 
progenitor has had time to become phase mixed predicts 
over-densities in places were the data show none. By con- 
trast, the dissolving, not-yet-phase-mixed scenario was 
able to account for the localization as well as reproducing 
the observed structure of the Aquarius stream. We there- 
fore suggest that the stream is dynamically young: the lo- 
calization could be explained as a recent disruption event 
of a progenitor whereby the stream has yet to occupy the 
available phase-space. The progenitor could either be a 
globular cluster or a dwarf galaxy, which may or may not 
have survived to the present day. We make no positive 
identification of with any globular clusters, though there 
could be a possible link with likely dwarf galaxy remnant, 
uj Cen, and the associated Kapteyn group. Follow-up 
high-resolution abundances would elucidate this possible 
connection. Further, more sophisticated simulations of 
Aquarius are required. This will enable a better under- 
standing of this interesting, new halo stream which places 
hierarchical formation right on our proverbial doorstep. 
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APPENDIX 

RULING OUT SAGITTARIUS DWARF DEBRIS 

We compared the kinematic properties of the Aquarius stream stars t o models of the Sa gittarius dwarf debris, using 
the prolate (q=0.9), sphe rical (q=1.0) and oblate (q=1.25) models of Law et al. (2005) (hereafter L05Fn We also 
used the triaxial model of Law et al. (2009) (hereafter L09) which uses a halo potential with c/a ~ 0.677o/a sa 0.83 
within 60kpcp^ We follow convention and use the parameter Aq (defined in Figure 1 of L05) in our plots, which is 
the longitude from the Sgr dSph in the plane of its orbit increasing away from the Galactic plane. Figure l9] plots Aq 
against Galactocentric radial velocity for the Aquarius stream stars and the different Law models. We highlight those 
stars from the models that have the following properties: 

• distance to sun < 15kpc 

• declination < 0° 

• Galactic latitude < 0° 

We are generous with the distance and Galactic latitude criteria to allow for uncertainties in the models (as well as 
the observational distances) . Figure |9] shows VgaX as a function of A© . Those stars that fit the above selection criteria 
and that have V^ai > — 50kms^^ are marked in green while those that have Vgai < — 50kms~^ are marked in blue. 
The reason for this delineation will become evident below. 

The velocities show a main feature at Vgai — +200 — 250kms~^, which corresponds to the leading arm, which 
the models predict is vertically streaming through the solar neighbourhood. This stream gets stronger going from 
the prolate towards the oblate models. This large signal is not pr esent in the RAVE da ta, as the heliocentric radial 
velocity is of the order of Vi 



Seabroke et al. 



(|2008 1 showed that there is no such 



-270kms~' or = -300kms 

large asymmetry detectable in the distribution of radial velocities for stars with / < —45". We do see, however, that a 
faint signal of stars is present for the oblate models with V^ai = — lOOkms"^, which is associated with an extra wrap 
of the leading arm passing south of the solar neighbourhood in this model. A value of V^ai = — 50kms~^ separates 
this extra, fainter wrap from the main leading arm component predicted by the oblate model. The triaxial model 
does not exhibit this feature and indeed the V^ai = +300kms~^ stream is weaker as the leading arm misses the solar 
neighbourhood. 

We concentrate on the oblate model with its possible curl of the leading arm fitting the Aquarius stream stars, as 
this is the only possible match. In Figure 10 we show the X-Z plane as well as the Vr-V^-Vz planes for the stream 
stars and the solar neighbourhood stars from the oblate model. The space velocities were calculated using the radial 
velocities and proper motions in the RAVE catalog as well as the distances derived via the two different methods 
(isochrone and RPMD). There is some overlap between the velocities from the two different distance derivations but 
in general we see that the space velocities are affected by the uncertainties in the distances. The RPM distances give a 
much tighter grouping in velocity for the stars and we take these results to be more indicative of the group's properties, 
plotting median error bars for these values. 

The first thing to note is that the V^ai < — 50kms~"'^ simulation particles do not fit the positions of the Aquarius 
stream stars in the X-Z plane. Even accounting for distance errors the distribution is strikingly different. Rather, the 
group stars tend to be aligned spatially with the Vgai > — 50kms~^ simulation particles. Secondly, while the Vr and 
Vz values for the Vgai < — 50kms~^ simulation particles are similar to that of the group, the values for very much 
differ from those of the Aquarius stream stars. For while the errors for the stream's values are larger than the other 
velocity components, a significant (~ 2<t) and systematic shift would be required in this component for the stream 
and V^ai < — 50kms^ simulation particles to agree. So while there does appear to be some overlap between this faint 
wrap and the group in a couple of variables, both the spatial distribution and the velocity distribution do not match. 
This is further borne out by the proper motions: the average value for the Aquarius stream stars is 28masyr~^ while 
for the Vgai < — 50kms~^ particles it is 4masyr^^. 
It is further wo rth noting that the ob late halo potential mo del docs not compare well to Sagittarius dwarf debris 

2006 ) data set traces dynamically old Sgr stream stars 
ark halos give different predictions. These data do not 



data. As noted by Fellhauer et al. (20061, the 
around the North Galactic C ap, where the ob 



Belokurov et al. 
ate a nd prolate 



favour the oblate model with Fellhauer et al. ( 2006 ) argui ng for a spherical dark halo, while L09 favo ur a triaxial halo. 
Furthermore, the absence of the Sgr stream iiear the Sun ( Seabroke et al.|2d08 Newberg et al.|2009 l is consistent with 



simulations of the disruption of Sgr in nearly spherical and prolate Galactic potentials. I'hus, the only model that has 
a passing resemblance to the Aquarius stream stars is the least likely of all those presented. On the strength of all the 
evidence, we conclude that the Aquarius and Sagittarius streams are unrelated. 



available from http://www.astro.virglma.edu/~srm4n/Sgr/ 



kindly provided before public release by D. Law 
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Figure 9. Longitude in the Sagittarius orbital plane, Aq, vs galactocentric radial velocity for the Aquarius stream stars (red points) 
compared to the Law models. Stars in the models with d < 15kpc, b < 0°, I < 0°, Vgai > — 50kms~^ are labelled green while 
d < 15kpc, b < 0°, I < 0°, Vgai < — 50kms^^ are blue. In this plot the blue points in the oblate model show a possible match to the 
Aquarius stream stars. 



The Aquarius Stream 
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Figure 10. XZ plane and V^, Vr and Vz plots showing the Aquarius stream stars and the oblate model of L05. The solid red points 
show the values for Aquarius stream stars using the RPM distances while the open red points use the isochrone derived distances. Median 
errors bars for the RPM distance derived values are also shown. The blue and green points are as in Figure [9] where we see that with the 
full 6D phase-space information the possible match with the blue points from the oblate model is ruled out. 



